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Abstract 

We construct the relationship between nonrenormalizable, effec- 
tive, time-reversal violating (TV) parity-conserving (PC) interactions 
of quarks and gauge bosons and various low-energy TVPC and TV 
parity-violating (PV) observables. Using effective field theory meth- 
ods, we delineate the scenarious under which experimental limits on 
permanent electric dipole moments (EDM's) of the electron, neutron, 
and neutral atoms as well as limits on TVPC observables provide the 
most stringent bounds on new TVPC interactions. Under scenarios 
in which parity invariance is restored at short distances, the one-loop 
EDM of elementary fermions generate the most severe constraints. 
The limits derived from the atomic EDM of 199 Hg are considerably 
weaker. When parity symmetry remains broken at short distances, 
direct TVPC search limits provide the least ambiguous bounds. The 
direct limits follow from TVPC interactions between two quarks. 



1 Introduction 



The search for physics beyond the Standard Model (SM) is a topic of consid- 
erable interest in high-energy particle physics. Concurrently, efforts are also 
underway to uncover signatures of new physics at low- and medium-energies 
using atomic and nuclear processes. In this respect, there exist tantalizing 
hints of new physics in results of neutron and superallowed nuclear /3-decays, 
which imply a value for \V u d\ differing from the SM unitarity requirement 
by two or more a [|TJ. Similarly, the weak charge of the cesium atom, Q w , 
measured in atomic parity-violation (APV) by the Boulder group, has been 
found to differ from the SM prediction by 2.5a (See, however, Ref. ||). If 
conventional many-body atomic and nuclear effects can be ruled out as the 
source of these deviations, the /3-decay and APV results imply the existence 
of new physics at the 1-10 TeV scale||, [|. This possibility has motivated a 
variety of additional atomic and nuclear new physics searches, including new 
measurements of the neutron /3-decay parameters, APV observables along a 
chain of isotopes, and parity-violating (PV) electron-electron and electron- 
proton scattering. 

One possible manifestation of new physics not probed by the aforemen- 
tioned experiments would be the existence of new low-energy interactions 
involving a single generation of fermions which violate time-reversal invari- 
ance (T) but conserve parity invariance (P). Such interactions are allowed 
in the SM when quarks of different generations participate. Recently, the 
first non-zero result for a AS = 1 T-violating, P-conserving (TVPC) observ- 
able has been reported by the CPLEAR Collaboration, which measured the 
K° — K° decay asymmetry || . The results are consistent with the value ex- 
pected from the measured CP-violating parameter e and the CPT theorem. 
No new physics is required to explain this result. In the AS = sector, 
a variety of direct searches for TVPC effects have been carried out. These 
efforts include studies of detailed balance in nuclear reactions^, 7-ray corre- 
lations in nuclear 7-decayj5|, five-fold correlations (FC) in the scattering of 
epithermal neutrons from aligned nuclear target s|10|, Oj , charge symmetry 



breaking (CSB) in np scattering [|12|, [13], [14|], and the J ■ (p e x p v ) in neu 



tron /3-decay P5| . Thus far, all studies have yielded null results. The limits 
from the purely hadronic reactions imply olt < few x 10~ 3 , where a? gives 
the ratio of TVPC nuclear matrix elements to those of the residual strong 
interaction. 
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Limits on AS = TVPC interactions involving light quarks may also 
be derived indirectly from results for atomic, neutron, and electron electric 
dipole moments (EDM's). As observed in Ref. [IB|, the presence of both a 



new TVPC interaction and a conventional PV interaction (e.g., in the Stan- 
dard Model) could conspire to generate a non-zero EDM, whose interaction 
with an external field violates both P and T. To the extent that PV radia- 
tive corrections to possible new TVPC interactions can be calculated, one 
can derive limits on new TVPC interactions from EDM results. Attempts to 
do so were first reported in Ref. 
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The calculation involved two external 
elementary fermions (e.g., two valence quarks in the neutron) and a one-loop 
Z-boson radiative correction to dimension seven, four-fermion, TVPC oper- 
ators. Two-loop effects, involving a single external fermion, for the electron 
EDM, d e , and neutron EDM, d n , were later studied in Ref. 117]. Naively, 



one might expect the most stringent bounds on new TVPC interactions to be 
derived at one-loop order from the experimental limit on the atomic electric 
dipole moment of mercury, d A ( 199 Hg), since the latter is nearly two orders of 
magnitude more severe than the bound on d n . It was argued in Ref. []I7], 
however, that the two-loop effects in d e and d n generate considerably more 
stringent bounds than do the results of d A . Subsequently, the authors of Ref. 



lq| recast the analysis of Refs. [0, [17]] into the framework of low-energy ef- 
fective field theory (EFT). It was argued in Ref. |18| that the results of Ref. 



IT] imply bounds on new TVPC interactions in excess of those presently 



achievable with direct TVPC searches by several orders of magnitude. These 
conclusions have had a discouraging effect on further direct TVPC searches. 
Recently, it was argued that the conclusions of Ref. |17] are inconsistent 



with the separation of scales underlying EFT PD| |. In brief, the argument is 
as follows |H 



Let A TVPC denote the mass scale below which use of an EFT 
involving nonrenormalizable TVPC operators makes sense. One may expand 
the EDM of an elementary particle, neutron, or atom as 



d = (3 5 C 5 
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where the Cd denote the set of a priori uknown coefficients of dimension d 
nonrenormalizable operators in the effective Lagrangian, the fld are calculable 
quantities arising from loops or many-body matrix elements, and M < A TVPC 
is a mass scale associated with the appropriate dynamical degree of freedom 
in the EFT. The Cd parameterize one's ignorance about the short-distance 
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(p> A TVPC ) dynamics of the new time-reversal violating physics. The first 
contributions from new TVPC interactions appear in the C7. 
One may now consider Eq. (|I]) under two scenarios: 

Scenario (A) Parity symmetry is restored at some scale fi < A TVPC . In 
this case, all of the coefficients C5 and Cg must vanish at tree-level in the 
EFT since parity invariance holds at short distances. Consequently, the first 
contributions to the EDM arise from loops involving the TVPC Cj operators. 
Since M/A 

tvpc < 1) these contributions presumably dominate the remaining 
terms in the series. Hence, one may use experimental EDM limits to constrain 
CV/A| VPC . As shown below, the limits obtained from EDM's under this 
scenario vastly exceed those obtainable from direct searches. 

Scenario (B) Parity symmetry is restored at fi > A TVPC . In this case, the 
C5 and C*6 do not, in general, vanish at tree-level in the EFT since both PV 
and TV interactions take place at short distance. There exists no reason 
to assume they fail to conspire in generating the lowest dimension TVPV 
effective interactions. Consequently, the TVPC interactions do not generate 
the leading contribution to the EDM as in Scenario (A). Without independent 
information on the C^^ one cannot use the EDM as a direct handle on the 
TVPC Cj terms. The latter may be more or less suppressed relative to the 
lower dimension contributions depending on the size of M/A TVPC . Since one 
has no a priori information on M/A TVPC , one can say very little about the 
importance of TVPC contributions. For the sake of argument, one might 
assume M/A TVPC << 1 so that the first term in Eq. (H) dominates. In this 
case, the low-energy effects of TVPC interactions would be negligible. In 
the more general situation, however, one would have to use direct TVPC 
searches to constrain the new TVPC interactions under this scenario. 

The analysis of Refs. |T7| , [U| implicitly assumes Scenario (A). The EDM 
calculations performed by these authors, however, do not display the proper 
scaling behavior which follows from EFT. The reasons for this failure 
are discussed in Ref. |21J and summarized below. It was also shown in Ref. 
[pCfl that there exist additional TVPC operators, not considered in Refs. 



bj, [L7|, |18|1 , which contribute to the elementary fermion EDM at one- loop 
order. Under Scenario (A), these one loop effects yield the most stringent 
constraints on the size of TVPC effects. 



In what follows, we extend the analysis of Ref. [2(| to include many-quark 
TVPC contributions to the neutron EDM - first studied in Ref. [UBI - and 
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to atomic EDM's. We concentrate on Scenario (A), since under Scenario (B) 
one cannot use EDM's to derive unambiguous information about TVPC new 
physics. In the case of the neutron EDM, we complete the one-loop analysis 
of Ref . jTjJ , including additional diagrams required by electromagnetic gauge 



invariance. We show that the impact of these new diagrams is as large as 



the one-loop effects considered previously in Ref. |Tj|. We also compute 



tree-level contributions arising from dimension seven TVPC operators not 



considered in Ref. 16 



In the case of atomic EDM's, we consider the situation in which they arise 
from purely hadronic TVPC interactions in the nucleus. Traditionally, the 
effects of non-leptonic T-violation in nuclear and atomic processes have been 
analyzed using collective degrees of freedom (mesons and baryons), rather 
than fundamental quark-quark or quark-gluon interactions. The T-violating 
effects are characterized by hadronic coupling constants which may be related 
to the underlying quark and gluon T-violating interactions using standard 
hadron structure techniques. In this context, two hadronic effects are of 
interest: (a) the presence of a purely TVPC meson- nucleon interaction and 
(b) the presence of a TVPV meson-nucleon interaction. 

The leading "long-range" TVPC effect arises from p-meson exchange, 
where the TVPC pNN vertex is characterized by a coupling strength g p and 
an interaction ]nj 

/™ = iy/2g p f p ^-N < x^(T-d xP +-T + d x p-)N , (2) 

where f p = 2.79 and Ky = 3.70. A time-reversal violating parity violat- 
ing (TVPV) nuclear effect arises when the second vertex in the exchange is 
parity-violating. Alternately, a TVPV atomic moment can be generated by 
a TVPC nuclear p-exchange and the PV exchange of a Z-boson between the 
nucleus and atomic electrons (see Fig. [1]). 

Similarly, the longest-range TVPV effects generally arise from 7r-exchange. 
In this case, the relevant TVPV irNN couplings are g$p', where the super- 
script denotes the isospin channel and correspond to the interactions 



nTVPV _ N 



^V-vf + ^V (3) 



+# / (3r z 7r -r-vr) N 
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Non-zero values of g^>' may arise from either new TVPC interactions plus 
weak radiative corrections, or from more conventional TVPV interactions, 
such as the #-term in the QCD Lagrangian. The latter have been considered 
extensively elsewhere [|H|, and we concentrate on the former. 



Upper bounds on \g p \ have been derived from a variety of T- violating, P- 
conserving experiments, including the studies of detailed balance in nuclear 



reactions^, neutron transmission through an aligned Ho target[10], and 



CSB terms in the np scattering cross section [[12]]. In addition, measurements 
of d n and d A yield bounds on g p , when PV and TVPC interactions conspire 
to generate an EDM. Under Scenario (A), the EDM limits on g p may be 



stronger than those obtained from TVPC nuclear processes |27]]. Regarding 
the g)p', experimental EDM limits yield the most stringent bounds. As 
discussed in detail in Ref. |L9|], d n and d A ( 199 Hg) are sensitive to different 



linear combinations of the g^p'. Roughly speaking, the EDM upper bounds 
on the \gjp'\ are of the order of a few xlO~ n . 

In this paper, we relate g p and the g(p' to the underlying TVPC quark 
and gluon operators, and use limits on the hadronic couplings to infer limits 
on these underlying interactions. We then compare these limits as well as 
those obtained from the one-loop many-quark contributions to d n with those 



obtained from the EDM's of elementary fermions as in Ref. [20|. We find: 



1. Under Scenario (A), the one-loop elementary fermion EDM studied in 
Ref. |2(J produces the most stringent limits on new TVPC interactions 



when applied to the electron and neutron. Many-quark effects d n or 
dA are suppressed by at least five orders of magnitude. 

By making certain naturalness assumptions, one may use experimental 
EDM limits to derive lower bounds on A TVPC . Under Scenario (A), one 
infers from d e (single-quark d n ) limits that A TVPC > 260 TeV (~ 110 
TeV) if the new TVPC physics is strong. The corresponding bounds 
derived from many-quark effects in d n and d A ( l99 Hg) are at roughly 
1000 times weaker. In order for d A ( 199 Hg) to provide competitive limits, 
the precision of the atomic EDM experiments would need to improve 
by roughly nine orders of magnitude^. 



3. Under Scenario (A), one expects ay < 10 



-15 



however, at some point, improved precision in d A would improve tighten the bounds 
on cL . 
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4. For Scenario (B), experimental limits on g p derived from FC in neutron 
transmission and CSB in np scattering produce the strongest limits on 
TVPC new physics. In terms of mass scales, these limits give A TVPC > 1 
GeV for new strong TVPC physics. An improvement of six orders of 
magnitude in experimental precision would bring this lower bound up 
to the weak scale. 



The analysis leading to these conclusions is presented in the remainder 
of the paper as follows. In section 2, we review the framework of EFT for 
new low-energy TVPC and TVPV interactions. In section 3 we illustrate 
the application of this framework by considering the EDM of an elementary 
fermion, as discussed in Ref. |2(J. In section 4, we discuss the renormal- 
ization of effective TVPV four-fermion operators arising from PV radiative 
corrections to effective TVPC interactions. Here, we take particular care to 
implement electromagnetic gauge invariance. The latter implies the existence 
of additional contributions to the EDM's of composite systems not considered 
in Ref. Hi~6| . In section 5, we relate the effective TVPC and TVPV oper- 
ators to g p and the gft', respectively, using the quark model, factorization, 
and current algebra. We also compute new many-quark contributions to d n 
generated by new operators - including those required by gauge invariance 



- not considered in Ref. |T(|. In section 6, we compare the implications of 



the many-quark contributions for the scale of new TVPC interactions with 
those obtained from the study of Ref. [20|. We also consider the limits on 



g p and the g^V obtained from atomic EDM limits and direct searches and 
the corresponding implications for new TVPC interactions under the two 
scenarios outlined above. Section 6 summarizes our conclusions. 



2 Effective field theory and new TVPC inter- 
actions 

Herczeg et al have shown that TVPC interactions between quarks cannnot 
arise from tree-level boson exchange in renormalizable gauge theories [§]. 
Such interactions could, however, be generated by higher-order or non-perturbative 
effects. Whatever new physics produces P-conserving T-violation among 
light quarks and gluons must be characterized by some heavy mass scale, 
A TVPC . Given that the underlying renormalizable gauge theory for P-conserving 
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T-violation is not known, it is natural to consider the low-energy conse- 
quences of such a theory in the context of an EFT valid below the scale 
A TVPC . Letting C NEW denote the effective, low-energy Lagrangian for new 



physics, we follow Ref. [18] and expand in inverse powers of A 



TVPC- 



C-new — £4 + — £5 + -r^ Cq + — £7 + • • • , (4) 

^TVPC ^TVPC ^TVPC 

where the subscripts denote the dimension of the operators appearing in each 
term and where 

4i = £C}0j (5) 
k 

with the sum running over k complete set of dimension d operators {O^}. 
The lowest dimension T- violating operators appear in £5. These operators 
are TVPV only. Of particular interest here are the electric dipole fermion- 
gauge boson interactions: 

OV = -\i>fW^ f F^ (6) 
0{ Z = -faa^f Z>" (7) 

0{ 9 = -^/<VT*AV/ G a , (8) 

where F^, Z^, and G a ^ denote the photon, Z-boson, and gluon field 
strength tensors, respectively, and the A a are the Gell-Mann matrices. 

The term £ 6 contains the lowest-order TVPV four-fermion operators, 
which include 

Ola = i^f^f^f'lBlpf (9) 

oW = i^ f X a ^ r X a i^ f (10) 

oW = ijfv^rfrv^f, (11) 

O^' = ^AV^V/^/'AV^tsV/' (12) 

and so forth. 

The lowest- dimension TVPC interactions arise in £7. Here, we consider 
the following three: 

O f 7 [' = z$/75<v(^ + ^)^/'7^75^'+ h - c - ( 13 ) 
Oft = ^Xa^fGiTK (14) 
Of? = $ f a^ f Z^F» . (15) 
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The operator 0{[ was first considered in Refs. [I7[], while O^ was 

introduced in Ref. The interaction Of c 7 was subsequently considered in 
Ref. im 



The d = 7 Lagrangian also contains several TVPV operators. Among 
those relevant to us are 

Off = ^CD^-^^Mfl^f (16) 

off = Mrfr-fisCD^^r (17) 

off = ii> n v i> f i>r°^(D» + ^Mf> ( 18 ) 

off = ijfYi^rtra^CDv-^M' (19) 

off = jffMr'ft&v-itM' (2°) 

Off = ijffCD" + ^)^/#^ 75 ^ (21) 

Off = i^ n v l50»-^Mfi>fi°^f> ■ ( 22 ) 

The specific forms of other d > 5 TVPV operators and d > 7 TVPC operators 
are not relevant to the present discussion, so we do not list them explicitly. 

A key ingredient underlying the expansion of Eq. (f|) is a separation 
of scales and an associated power-counting scheme. Specifically, the contri- 
bution to a T-violating observable from any physics associated with scales 
\i > A TVPC is contained in the operator coefficients, These short-distance 
contributions are not calculable since the underlying theory (e.g., renormaliz- 
able gauge theory) responsible for TVPC effects is not knownF]. Consequently, 
the Cd can only be determined from experiment. Contributions from physics 
having /j, < A TVPC live in loops and many-body (e.g., hadronic) matrix ele- 
ments containing the non-renormalizable operators Od and physical degrees 
of freedom having masses and momenta less than the scale A TVPC . Only 
these "long-distance" contributions can be computed using the EFT. 

As a consequence of this scale separation, one obtains a systematic power- 
counting scheme by which to organize contributions to any low-energy T- 
violating observable. If d is the lowest dimension of an effective operator 
which contributes to such an observable A T ~ ODD , then 

jt-odd _ Cd (_L\ d ~\... j (23) 



2 If this underlying theory were known, one would have no need to employ the expansion 
in Eq. (Q) in the first place 
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where p denotes a mass or momentum smaller than A TVPC . To the extent 
that p « A TV pc, contributions from higher- dimension operators will be 
suppressed relative to those from Od by powers of (p/A TVPC ). In general, one 
may thus truncate the expansion of J?-° DD at the first or first few orders 
in (p/A TVPC ). In our analysis of the EDM, we find p is one of the following: 
elementary fermion mass, mf, weak gauge boson mass, M z ; QCD scale, 
A QCD ; inverse hadron size, l/r HAD ; and long wavelength photon momentum, 

Any renormalization of the Od by loops involving any one of the Od and, 
e.g., degrees of freedom appearing in £4 must be carried out in a manner 
which preserves the EFT scale separation. In particular, only intermedi- 
ate states with momenta and energies below A TVPC must contribute to loop 
integrals. Higher-momentum (p> A TVPC ) states have been effectively inte- 
grated out in arriving at the expansion in Eq. (|j). Consequently, we regulate 
all loop integrals using dimensional regularization (DR), which preserves the 
separation of scales. Since the subtraction scale \x arising in DR only appears 
logarithmically in any regulated amplitude and never as a power, the use of 
DR does not alter the EFT power counting described above. We emphasize 



that the cut-off regulator used in Refs. [16, 17] does not preserve the EFT 



scale separation. In those analyses, loop integrals were cut off at momenta 
p ~ A TVPC . Consequently, in Refs. |16|, [T^] the renormalization of the EDM 
5a due to loops involving any d > 5 operator scales as 1/A TVPC and not as 
(P/A tvpc)^ 5 x (1/Atvpc) as implied by EFT power counting. As we show 
below, this loss of power counting prevents one from deriving any information 
about the d = 7 operators from experimental EDM limits. In this respect, 



our conclusions differ dramatically from those of Refs. [16, [17|]. We illustrate 
this point in the following section. 



3 Elementary fermion EDM 

The application of EFT for new TVPC interactions to the EDM of an ele- 
mentary fermion was considered in Ref . |2(J . In what follows, we summarize 



the arguments of that analysis, as they illustrate the general principles of 
EFT to be used in the remainder of this study. To that end, we first observe 
that if, as in Scenario B, both new TVPC interactions and PV interactions 
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(e.g., in the Standard Model) exist at momentum scales p^A TVPC , then there 
exists no reason to assume that the coefficients of the TVPV effective opera- 
tors in Eq. vanish at tree level. Although we have no detailed knowledge 
of the dynamics of short- distance TVPC and TCPV interactions, nothing 
prevents their conspiring to generate non- vanishing low-energy TVPV inter- 
actions. In particular, the coefficients of the d = 5 electric dipole operators, 
C 5 /7 should be non-vanishing at tree-level unless some fortuitous fine-tuning 
of the short-distance TVPC and TCPV interactions occurs. 

The situation here is analogous to the chiral expansion of the octet baryon 
magnetic moments. In the latter case, the leading order contribution occurs 



at tree level from the d = 5 magnetic moment Lagrangian [ 22j : 

e 
2A 



{b+Tr (B V S?{\ 3 + X s /y/3, B v }) (24) 



+6_Tr(^[A 3 + A 8 /v / 3,£q) 



where the B v are the octet baryon fields for states of velocity v a , S a is the 
spin operator, and A x is the scale of chiral symmetry breaking. The tree 
level relationship between a baryon magnetic moment and the low-energy 
constants b± is 

, (25) 

where m B is the mass of the baryon, "a" denotes its SU(3) indices, and b a is 
the appropriate combination of the b±. Since the baryon magnetic moments 
are typically of order unity, the tree-level relation implies that the b± are also 
of order unity. Alternatively, one may use Eq. ( PS| ) to estimate the scale A x . 
If the low-energy constants in Eq. ([24]) are of order unity, the one must have 
A x ~ m B ~ 1 GeV. Since the chiral symmetry of pionic interactions implies 
that A x = 4ttFk ~ 1 GeV, the tree-level magnetic moment relation produces 
a self-consistent value for the scale of chiral symmetry breaking when the 
leading low-energy constants are chosen to be of order one. 

In a similar way, one may use the tree-level relation between the EDM 
and the coefficients of the appropriate d = 5 operators to estimate the scale 
A TVPC . This relation is 

d, = -p- • (26) 

J-'-TVPC 
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We follow a standard convention for parameterizing the strength of new 
physics interactions and take C5 7 = AnK 2 e. Using the present limit for the 
EDM of the electron \d e \ < 4 x 10~ 27 e-cm[|I], 0, one obtains from Eq. (fZ5|) 
the limit A TVPC > 10 14 k, 2 GeV. Thus, if the new TVPC physics is "strong" 
(k, 2 ~ 1), one obtains a tremendously large value for the corresponding mass 
scale. 

For both EFT's in Eqs. (|], [24]), loop corrections involving light, dynami- 
cal degrees of freedom modify the tree- level relations in Eqs. (^,^). Single 
pion loop amplitudes, for example, generate corrections to isovector mag- 
netic moments of 0(m n /A x ) relative to the tree-level contribution. The ir 
loops are quadratically divergent, yet generate a finite contribution to the 
d = 5 magnetic moment operator. Power- counting implies the appearance of 
one additional mass factor in this finite contribution. When DR is used to 
regulate the integral, this mass factor becomes m n , resulting in the m n /A x 
suppression relative to the tree level relation. This scaling behavior of the 
loop contributions, which follows from the EFT separation of scales, pro- 
vides for a chiral expansion in powers of p/A x (where p < A x ) which may 
be reasonably truncated at any order. Similarly, if the EFT of Eq. (£|) is 
well-behaved, one would expect the loop corrections to the relation in Eq. 
( |26|) to be suppressed by powers of p/A TVPC , where p < A TVPC is a mass scale 
associated with the dynamical degrees of freedom in the C NEW . 

For purposes of illustration, we assume that the only dynamical degrees of 
freedom operative below A TVPC are those appearing in the Standard Model. 
Under Scenario A, additional degrees of freedom - such as the right-handed 
neutral gauge boson in left-right symmetric theories - would also generate 
loop contributions. These additional degrees of freedom would be required 
in order for parity symmetry to be restored for fi < A TVPC . As noted below, 
however, the study of only the standard model contributions yields conser- 
vative upper bounds on the d = 7 TVPC operators. In order to avoide 
introducing model-dependence associated with parity restoration scenarios, 
we restrict our attention to these SM effects. In this case, the leading loop 
corrections are generated by PV Standard Model radiative corrections to the 
d = 7 operators 0^ a - c . These corrections have been computed in Ref. p(J for 
O-j^c using the diagrams in Figs. || and §. All of the loop amplitudes corre- 
sponding to Figs. |2] and |] are superficially quadratically divergent. The loops 
are regulated using DR and the poles removed by the corresponding coun- 
terterm in the renormalized C[ in the MS subtraction scheme. The closed 
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fermion loop subgraph in Fig. is nominally linearly divergent and corre- 
sponds to the Adler-Bell-Jackiw anomaly diagram. In this case, the vector 
current insertions arise from the 7 and Z couplings to the internal fermion, 
while the axial vector insertion arises from Oj a . Since the EDM operator is 



linear in the photon momentum q^, we follow Ref. |L7| arid retain only the 
terms linear in arising from this sub-graph. Denoting its amplitude T AtA ", 
we choose the loop momentum routing to satisfy q^T^xa = = k x T pXa , where 
q p and k\ are the photon and Z-boson momenta, respectively. The result is 
the usual anomalous term in (q + k^T^xa- To linear order in q, there exist 
three structures which satisfy these vector current conservation conditions: 

A^ Xa = k ■ qk p e^ xpQ - k»e aXpa q a k p (27) 
B nXa = k 2 q„e Xf * va - k x e p ^ a k p q v (28) 
C» Xa = k a e apxp q a k p . (29) 

The loop integrals for T^ Xa are nominally linearly divergent. As a check 
on the calculation, we regulate the integrals using two different regulators - 
Pauli Villars and DR - and obtain identical results in each case. The result 
is finite: 

r*. = + 5B ^ _ ac^] I ix .(so, 

It is straightforward to verify that when the corresponding term linear in k 
and quadratic in q, the photon momentum^], is added to the expression in 
Eq. (|30f) and the divergence (q + k) a T p x a computed, one obtains the finite, 
textbook result for k 2 = q 2 = |23|]. 

The closed fermion loop of Fig. |^b contains axial vector insertions from 
Oj a and from the coupling of the Z-boson to the internal fermion. The 
external fermion line contains the vector current Z-fermion coupling. This 
sub-graph diverges quadratically and must be renormalized by the appro- 
priate MS subtraction before the remaining loop integration is performed. 
In all cases, the amplitudes are infrared-finite. Consequently, we follow Ref. 
[|P7|j and neglect the fermion mass dependence entering the loops. 

At leading-log order, the results are 

C/~eC-(^) 2 (^)4-A,)lnM (31 ) 
\A TVPC / KswCwS Vlo7r/ ^ 



'Obtained from Eq. (BOJ) by k <-> q, /i «-» A 
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for the one-loop contribution in Fig. |2|, which contains 7c ||24j| . The loops 
in Fig. |3|, which contain the four-fermion operator 7a , yield 




for the amplitude of Fig. |3|b. Here, sw and cw denote the sine and cosine 
of the Weinberg angle, respectively, g{r and g{ denote the vector and axial 
vector couplings of the Z-boson to fermion /, and Qf is the corresponding 
fermion electromagnetic (EM) chargeQ 

As expected from general considerations, the results in Eqs. (|3~Tl-[32l) dis- 
play the (p/A TVPC ) 2 suppression relative to the tree-level relation in Eq. (|26|), 
where in this case p = M z . Thus, if the C 7 have natural size 

C 7a = 4™ 2 (34) 
Airc\ 

C 7c = AnK 2 eg= C 7a (35) 

sw 

and if A TVPC is determined from the experimental limit on d e via Eq. (|26|) , 
then the loop corrections to the tree-level EDM will be suppressed by more 
than 25 orders of magnitude for new strong interactions (/t 2 ~ 1). In this 
case, tree-level dominance of the EDM implies that EDM limits do not pro- 
vide direct bounds on the d = 7 operators appearing in Eq. (||). The situation 
here is analogous to that of the baryon magnetic moments. Since the latter 
are generally dominated by the tree-level term, they cannot be used to de- 
termine the meson-baryon couplings which enter the one-loop, sub-leading 
contributions. Instead, the meson-baryon interaction must be determined di- 
rectly from, e.g., tcN scattering and the results used as input into the chiral 
loop corrections. 

The experimental EDM limits could be used to constrain the d = 7 
TVPC operators directly if - as in Scenario (A) - C*/ vanishes at tree-level 

4 We have not considered loop effects involving 0!j% . Since they contribute to the EDM 
at two-loop order, however, we expect them to be no more important than the two-loop 
contributions involving 7a . 
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in the EFT. Such a situation could arise under scenarios, such as left-right 
symmetric gauge theories, in which parity symmetry is restored at some 
scale well below A TVPC . In this case, there would exist no short- distance PV 
interactions to conspire with the new TVPC physics in generating a tree-level 
C( . The leading contribution to an elementary fermion EDM would then be 
given by the results in Eqs. (|31"1 - |3~3"D , with M z replaced by, e.g., the scale 
of parity-restoration (such as the mass of a right-handed gauge boson), and 
with the appropriate combinations of couplings. Since the low-energy scale 
(e.g., M z ) enters quadratically, a conservative upper bound on the C-jj IS? TVPC 
for this scenario can be obtained by using the Standard Model results given 
above. The most stringent constraint results from the one-loop amplitude 
in Eq. ( |5T] ) applied to the electron EDM. Using the parameterization of Eq. 
(p5|), we obtain A TVPC > 260 k 2 ^ TeV. The corresponding two-loop limits 
are about a factor of five ~ (167T 2 ) 1 / 3 weaker. The one-loop constraints from 
the neutron EDM are also slightly weaker, given the somewhat less stringent 
experimental limits on d n p9| . 

The results of the foregoing analysis have important implications for low- 
energy direct TVPC searches in light quark systems. These implications are 
most transparent under Scenario A. In this case, EDM limits constrain the 
ratio Cij ' K\ VPC via the one- and two-loop results of Eqs. (|3T] - p3i ) . As noted 
in Ref. fll8| , one expects the ratio «t to scale as 

a T ~ C 7 (p/A TVPC ) 3 , (36) 

where p is a momentum characteristic of low-energy hadronic interactions. 
Taking p = 1 GeV/c and using the one-loop electron EDM results, one 
obtains a limit of olt ^ 10~ 15 . By comparison, the present direct TVPC 
search limits are considerably weaker: ax ^ 10 -3 . In short, under the parity- 
restoration scenario (A), the EDM results provide the most stringent bounds 
by many orders of magnitude. 

For Scenario (B) in which C( does not vanish at tree-level (e.g., PV 
persists at short distances), the situation is more subtle. In this case the EDM 
results do not provide direct constraints on the d = 7 operators. Nevertheless, 
one might argue from the EDM limits that the low-energy effects of d — 7 
operators should be considerably smaller than the present sensitivity of direct 



TVPC searches. Comparing Eqs. (31-33) to Eq. (36), we infer that the low- 



energy TVPC effects of the d = 7 operators should be suppressed relative to 
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the corresponding contributions to an EDM by 



) < CI /loop fax-tors) . (37) 



If, in addition, the d = 7 loop contributions are already suppressed relative 
to the tree-level EDM by many orders of magnitude, one would conclude the 
that corresponding effects in low-energy (pM GeV/c) TVPC processes would 
be even smaller - certainly well below the present direct search sensitivity 

Nevertheless, this line of reasoning is not airtight. If A TVPC ~ M z , it is 
conceivable that the tree-level and loop contributions to the EDM can be 
comparable in magnitude and that, due to possible cancellations, the mag- 
nitude of either term can be considerably larger than the EDM limit itself. 
An analogous situation arises, for example, in the chiral expansion of the 
isoscalar nucleon magnetic moment. In this case, the leading corrections to 
the tree-level relation in Eq. fl25p arise from kaon loops. Since and A x 
do not differ appreciably, the loop corrections are consirably larger than the 
isoscalar magnetic moment. A similarly large short- distance (tree- level) con- 
tribution is needed to cancel the loop effect and obtain the small isoscalar 
magnetic moment. In this case, the use of the isoscalar magnetic moment, 
together with power counting and naive dimensional analysis, to infer either 
the size of A x via the tree-level relation Eq. ( p5|) or the size of the one-loop 
effects would lead to erroneous conclusions. An independent determination 
of the strength of the loop contribution (e.g., of the kaon-baryon interac- 
tion) is needed. Should a similar situation obtain for the EDM, then direct 
TVPC searches would still be needed to ascertain the scale of the d = 7 
contributions. 

Before concluding our discussion of the elementary fermion EDM, we 
emphasize the differences between our analysis and that of Ref. Jl7|. In that 



work, a calculation of the amplitudes in Fig. 3 was used to try and estimate 
the size of the short-distance contributions. This estimate was implemented 
by regulating the two-loop integrals corresponding to Fig. [3] with a form 
factor of the type 

F (P 2 ) = (P 2 /A 2 TVPC ~ I)" 1 ■ (38) 

The use of this regulator causes the loop integrals to be dominated by con- 
tributions from intermediate states having p ~ A TVPC , thereby blurring the 
separation of scales crucial to the EFT expansion of Eq. (ffl). Consequently, 
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the two- loop results of Ref. [I7j scale, incorrectly, as 1/A TVPC rather than 
as 1/K% VPC . The corresponding implications for low-energy direct TVPC 
observables are, therefore, erroneous. 

[20 



More generally, as noted in Ref. 



the use of a cut-off regulator destroys 
the power-counting which justifies truncation of the EDM analysis at d = 7. 
This loss of a systematic expansion can be seen by considering the tower of 
operators 

0% 2n = ^ f D^ f (d 2 r^ fn ^ fl , (39) 

where n 



0,1, 



The two-loop calculation of Ref. 



]17fl may be repeated 

regulate the divergences, one may, as the calculation of Ref. JL7|, regulate 
the integrals with a form factor 



f ft 

by replacing the insertion of 7a by each of the operators in Eq. 



F n {p 2 ) = (p 2 /A 



2 

TVPC 



1) 



-(1+n) 



(40) 



The corresponding loop integrals will be the same as those in Ref. 11 
with additional factors of 



2 

TVPC 



(V/A 
p 2 /A 

1 + n(p 2 /A 



(p 2 /A 







2 

TVPC 



2 \n-l 
TVPC) 



2 

TVPC 

(p 2 /A 
(p 2 /A 



1 



2 

TVPC 
2 

TVPC 



l) n + n{p 2 /A 2 TVPC ) 

\r n — 



n-l 



but 

(41) 

(42) 
(43) 



appearing in the integrand. The first term (= 1) on the RHS of Eq. ([43]) will 
yield the same leading-log contribution as obtained in the calculation of Ref. 
|j~7ij . The remaining terms will generate sub- leading contributions, finite as 
A TV p C — > oo. Thus, at leading-log order, each operator in the tower will 

f f 1 

generate the same contribution, apart from the operator coefficient C 7 + 2 m so 
that the EDM will be proportional to 



7+2n 



(44) 



n=0 



In this case, there exists no reason to isolate the effects of the d = 7 operators 
from those of any other operator in the tower. All contribute with equal 
weight. It would be erroneous, therefore, to truncate the series at d = 7, 
as was done in Ref. 0], and to argue that the EDM limits constrain the 



magnitude of only one term in this infinite series. 
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As this example illustrates, the preservation of the scale separation is 
crucial to the power counting arguments which justify truncation of the 
expansion of the EDM at a given order. When DR and MS subtraction 
is used, for example, the contributions of each operator in the tower (j39|) 
will be supressed by successive powers of (M Z /A TVPC ) 2 . To the extent that 
M Z /A TVPC < 1, truncation at d = 7 makes for a reasonable approxima- 
tion. In the remainder of this study, we therefore to work with DR and MS 
subtraction in treating loop effects. 



4 Four-quark TVPV Operators 

The previous discussion considered the EDM of an elementary fermion. For 
a composite system such as a neutron, for example, one must also consider 
many-body contributions involving more than one quark degree of freedom. 
A generic contribution of this type is shown in Fig. |]. The operators which 
describe these many-quark effects include the TVPV d = 6,7 operators listed 
in Eqs. (|| fL6|) . As in the case of the sing le fermion EDM C{ 7 , the d = 6, 7 
operators may exist at tree-level in the EFT if parity is violated at sufficiently 
high scales (Scenario B). Similarly, these operators may be renormalized by 
PV radiative corrections to the d = 7 TVPC operators. As with O5 7 , the 
Oq ^ pv will be dominated by these loop effects if parity symmetry is restored 
for fi < A TVPC (Scenario A). In what follows, we compute the relevant loop 
effects. 

The leading corrections to the 0^ PV are generated by the set of graphs 

illustrated in Fig. [5], where the operator inserted is 0{[ . The diagrams of 
Fig. |5|a were considered previously in the study of Ref . jTj| . The diagrams 
for Fig. [5|b, which are required by electromagnetic (EM) gauge invariance, 
were not included in that analysis. The inclusion of these graphs is needed in 
order to obtain the pieces of the d = 7 TVPV operators containing the pho- 
ton field. As we discuss in Section 6, the contributions from these 7-insertion 
diagrams to the neutron EDM are numerically as large as the contributions 
arising from the graphs of Fig. |5]a. The reason for this situation is relatively 
straightforward. The diagrams in Fig. |5]a yield the pieces of the d = 7 oper- 
ators containing a derivative, as well as contributions to the d — 6 operators 
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proportional to one power of quark massf]. When the derivative operator acts 
on quarks inside the hadron, the result is of order <A QCD . The resulting 
contribution to the neutron EDM, as in Fig. |]a for example, will therefore go 
as A QCD / AM , where AM is the mass difference between the neutron and one 
of its excited states (e.g., an unbound ir~p pair). The 7-insertion diagrams 
of Fig. ||b do not produce such derivative operators, and their correspond- 
ing neutron EDM contributions contain no A QCD /AM factors. To the extent 
that Aq CD I AM is of order one, the magnitude of the two sets of contributions 
should be comparable. 

The amplitudes, -M|, of Fig. |5] are logarithmically divergent. As before, 
we regulate the loops using DR and define the finite results using MS sub- 
traction. At leading-log order, the sum of all twenty six diagrams yields the 
following linear combination of the 0^ PV : 



M " = £3^b; log (w){ 6m '' 9 ^' a 



/ J'/nff 



(9a + 9a 



+9v- 



2 



30£f + o{'/ 



+m fg y v 0L f + (g{g f v ~ 1Qa9v ~ ^Wv) O f 7 '/ 



gWvOti + lgWvOf/} ■ (45) 



The expression in Eq. ([45]) is obtained by keeping the external fermion lines 
off shell. Doing so allows us to identify uniquely the contributions of the 
d = 7 TVPV derivative operators and verify the gauge invariance of the 
overall result. 

In order to compare our result with the calculation of Ref . |TB| , we use the 



equations of motion and let the quarks go on shell and convert to momentum 
space. We obtain 

CP' a / /x 2 \ r f f - 

M s = A 3 a o 97r ,2 r 2 lo s( T7T ) [terrif'girg'AUflsUfUf'Uf' 

-2im f ,g f v (g s A + U f U f U f j 5 U f 



5 Since the d = 6 operators do not preserve chirality, they cannot be induced dynami- 
cally by massless quarks. 
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+im fl g s ^g s v U f^a^U f U ftJ^U, f 
+i Wa9v ~ 9a9v + ~9a9v) Ufjsip'f + PfYUfUf^Uf 



+ MU f a^U f U f {p' fl +p r ) u7fll5 U f 



(46) 



where Uf = U(j)f) and Uf = U (p'j) are the spinors for incoming and outgoing 
fermions /, respectively. 

The first three of the terms on the RHS of Eq. (|46D are identical to 
those appearing in ||16|| . Our coefficient of the fourth term, however, differs 



from the corresponding expression in ||16|| , and the fifth term does not appear 
in that work at all. We trace part of the difference on the fourth term to 
diagrams in Fig. 5a where the Z° boson connects to initial and final quarks 



of the same species. It is unclear from the discussion of Ref. |16| whether 
those authors included this class of diagrams. Given that our sum of the 
amplitudes for Figs. 5a and 5b satisfies a gauge invariance self-consistency 
check, we are confident in our result. 



5 Hadronic Matrix Elements 

The operators obtained in the previous sections can be used to compute 
contributions to the EDM of the neutron and of neutral atoms. Doing so 
requires that one calculate various hadronic matrix elements of two- and 
four-quark operators. A first principles treatment of these matrix elements 
in QCD goes beyond the scope of the present study. Moreover, since we seek 
only to derive order of magnitude constraints on new TVPC interactions and 
not to obtain definitive numerical results, it suffices to draw upon various 



approximation methods. To that end, we turn to the quark model [p5 
factorization, and chiral symmetry. 

Below, we estimate a number of different matrix elements relevant to the 
neutron EDM and couplings g p and g^: 

1. The contribution to d n from quark EDMs (Fig. 6) 

2. The relationship between g p and d = 7 TVPC operators (Fig. 7a) and 
the relationship between g p and d n (Fig. 7b). 
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3. The contribution to d n from the four-quark/photon TVPV operators 
appearing in Eq. (f45| ) (Fig. 4b). 

4. The contribution to the from the purely hadronic terms in Eq. ( fj5|) 
(Fig. 8a) and the relationship between the and d n (Fig. 8b). 

5. The tree-level contribution from 0^2 (Fig. 9). 



A. Two-quark TVPV Operators 

Relating the EDM of a constituent quark to that of the neutron using 
the quark model is a straightforward exercise. As shown in Ref. |[20|| , this 
relationship is given by 



A 



TVPC 



I d 3 x (V + \e) , (47) 



where u and I are the upper and lower component quark model radial wave 
functions, respectively. The integral in Eq. (|47j ) can be estimated using the 
wave function normalization condition 



d 6 x [u 2 + t)=l (48) 
and expression for the axial vector charge 

g A = I [ d'x U - \e] , (49) 



3 J V 3 

where g A ~ 1.26. From Eqs. p3-H9|) one obtains 



B. Four-Quark TVPC Operators 

Deriving the relationship between the effective hadronic coupling g p and 
the four-quark d = 7 TVPC operators (Fig. 7a) requires more thought 
than in the case of evaluating two-quark matrix elements. For simplicity, we 
focus on the four-quark operator Oj[ . We make a simple estimate using 
factorization. Doing so requires use of the Fierz transformed version of this 
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operator, since the interaction in Eq. (0) involves only p . The Fierzed form 
of the operator is 

°7a = ~{^f 9 v^>fi>rlv^S + ^ ftp f'^ flu 9 V t/J f 

3-4— - 3- > 

+|^/7^ d vipyipftpf + -^ipn^fipf d v ip f 

1 _ ^ _ % — — ¥ 

+-ipf d ujfj.ipfipf^^ipf + ~£pfl pip pip f>°nv 9 v if) f 

1 _ ^ _ % — — > 

—gl>f d v^vipf'ipf'Juipf - -ipfa^f'tPf'ln 9 u^f 

3-4— - 3- - — > 

+-ipf d vlsipf'ipf'lvl^f + -ipfizipf'ip'fiviz d v ipf 

3-4— - 3- - — > 

-ji>f d uiui^f^Pfi^f - -^Pfivi^r^rivih d „ipf 

+ -^f 9 vl5^nu1pf'^f'l5ln1pf + -fPflSVvulpfilpf'lSlii'd V0/ 

1 _ 4 _ % — — > 

+jipf d u^nipfipf^a^ipf + -if} f m '^ip Tip /'75<V 9 „ip f 

(51) 

In the factorization approximation, one makes the replacement 

(N'\q t 1 O l q ] 2 q 3 2 2 q l 1 \Np) -> (0\q l 1 O 1 q J 2 \p}(N'\q J 2 O 2 q l 1 \N) , (52) 

where N and N' denote nucleons, is the field for a quark of flavor a and 
color i, and q^Oiq^q^O^ 1 denotes any of the products of quark bilinears 
appearing in Eq. (0). Note that since hadrons are color singlets, one has 

(Ol^Oxg/lp) = \8iM& X °rth\p) ( 53 ) 

(N'\q 2 2qi l \N) = l -5 i3 (N'\q m 2 2 q m l \N) , (54) 

where repeated indices are summed over. Hence, each factorization contri- 
bution contains a factor of 

-5a x -5n = - . (55) 
3 v 3 ji 3 v } 

Since any pieces of the interaction above which involve the 75 will not give 
rise to a p meson- vacuum matrix element, we are concerned with only the 
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first eight terms of the Fierzed interaction. We may also reduce the number 
of terms to be evaluated using the equations of motion. Of the resulting op- 
erators, we keep only those containing no powers of the quark mass, since the 
latter generate significant suppression factors. In the case where / represents 
an up quark and /' represents a down quark, the remaining structures are 

1 «_ i > i <— - i ... - — > 

-u d v dd^ u u + -u-y u dd d u u — -ud u -y p dda^u + -ua^ v ddj^ d v u. (56) 

In the factorization approximation, the matrix elements required are 

(0\ut^d\p')(n\dYu\p), 

(0\u^d\p-)(n\d3^u\p) (57) 

(0\ut^d\p-)(n\da^u\p), 

(0\w^d\p-)(n\d 1(i t u \p), 

plus the corresponding matrix elements for p + n —>■ p (note that in Eq. (|57|) 
the color indices have been suppressed for simplicity). A detailed evaluation 
of these matrix elements appears in Appendix A. The resulting TVPC pNN 
Lagrangian is 

c = ^^W^_i|E) r m 

+^ (°- 176 ^ " °- 122 ) <b>°*\ (r-pt ~ r + p-)N (59) 

It is customary to write the standard rho-nucleon Lagrangian and the 
TVPC rho-nucleon Lagrangian in the respective forms (see, e.g. Eq. @, 

C = V2f p N(Y + ^^^(r-pj + r+p-)N (60) 

C TVPC = V2g p f p ^-N^ x q x (r- P ; - r+ p-)N (61) 

where Ky is the anamolous isovector magnetic moment, m n is the mass of 
the nucleon and f p is the pNN coupling constant. By redefining the phase 
of the rho meson, we can eliminate the Dirac structure 7 M in Eq. (|58D . We 
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begin by writing Eqs. 



C 



pNN 



= V2f p N 
+V2f p N 



and 



0|) together as 



y (1 + 9) + z-^(l + iS^)o^g A 



'1 - 



where 



z'A, 



2m n 

' 

2m„ 



# mjpl.05 



6 Aj.y pc . fp 



fpK 

/> 
1.293 



r-p + pN 



B = - 



i a 



ud 
7a 



m. 



6 Aj,yp C fp 



0.176^ 

Rn 



Rn 

-0.122 



r + p-N (62) 



(63) 



(64) 



We then observe that since 1+iA/ f p ~ exp(iA/ f p ) and l+iB2m n / (f p ky) r 
exp[zJB/(/ p /cy)], the phases of the rho mesons can be redefined as, p+ — > 
p+ exp(— iAj f p ) and p~ — > p~ exp(iA/ f p ). This allows us to rewrite the total 
Lagrangian as 



p 

B2m„ 



2m. n 



fp^V 



A" 
fp, 



2m r 



T-ft + T + p-)N 
«""<ly{T-fi ~ r + p-)N 



The quantity g p is then 

2 ^ ' 



9p 



riud 
°7a 



m„Tn, 



A 3 
J^TVPC 



) fpKv 



1.05 _ 1.293 \ k v 
1T P ~~R^) 2m~ 



where 



1.05 1.293\ k v 



Rr, 



Rr, 



2m 



0.176 : 



N 



, Rn 

Rn 



N 



0.122 



0.176— 

Rn 



0.122 



-0.023 



(65) 



(66) 



(67) 



As calculated in Ref. ||27|| , g p may contribute to the neutron EDM via 
the loop diagrams of Fig. 7(b) as well as to the atomic EDM via process like 
those shown in Figs. l(a,b). The d n calculation of Ref. [27| included the 
introduction of form factors at the hadronic vertices in order to render the 
loop integrals finite. The result is 



dn h^NNgp-K-yfpQp^V £,/ « \ 



16^/2 



7rm, 



(68) 
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where g p7r7 = 0.4 and the PV pion-nucleon coupling h^NN is constrained by 
the PV 7-decay of 18 F to lie in the range: h n NN = (0.73 ± 2.3)g 7T , where 

= 3.8 x 10~ 8 characterizes the strength of the charged current AS = 
hadronic interaction. The function F depends on the masses appearing in 
the loop integral as well as the form factor cut-off parameter, A. 

We note that the use of Eq. (|68|) and the experimental limits on d n 
to derive bounds on TVPC interactions entails several ambiguities. First, 
the value of H-xnn measured in nuclei such as 18 F may differ from the value 
appropriate to the single nucleon or few-nucleon systems. Many-body nuclear 
effects may renormalize the long range PV NN interaction in such a way as to 
shift the value of the effective PV ttNN coupling from the value appropriate 
for Eq. (^). This ambiguity may be resolved by future experiments, such as 
the measurement of the n + p — > d + 7 asymmetry planned at LANSCE[E| 



Second, the use of a form factor to render the loop integral finite can 
introduce considerable ambiguity. The form factor chosen in Ref. [E7| was 
taken from the Bonn potential, with A = 1.4 GeV. One may just as well have 
chosen a cut-off given by the inverse size of the hadron. The variations due to 
this spread of choices can be significant [22]. Moreover, as argued in Section 



3, the use of cut-off regulators can render one's effective theory devoid of any 
systematic power counting, leaving it poorly defined. For these reasons, we 
will not use Eq. to derive limits on the d = 7 TVPC operators. 



C. Four-Quark TVPV Operators 

We consider the relationship between the d = 6 and d = 7 TVPV op- 
erators and the neutron EDM. We specify here to Scenario (A), in which 
case the TVPV operators arise entirely from PV radiative corrections to the 
TVPC operators, as in Fig. 5. First, we estimate the contribution from the 
pieces of the d = 7 operators containing the photon field. These contribu- 
tions can be understood diagrammatically as shown in Fig 4(b). Starting 
from expression in Eq. (|45|), we see that the TVPV 7- four-quark interaction 
can be written as 



C 



TVPV,EM 
eff 



eJ* ff A x 

^7 a 



a 



+ 



TVPC 
/' / 

9v9a 



32ns 



w c w 



log 



9v(9a+9a) ^fOx^rtf'lsl^r 
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+g(g{3i<ipf~/5il)f'ipf>j\ipf> 



A x . (69) 



For simplicity, we consider the case where / is an up quark and /' denotes 
a down quark. As in our estimate of g p , we can find neutron matrix elements 
using the quark model. We identify the appropriate quark model expression 
for the electric dipole moment, 

i(n\'\ J d 3 x 3 J eff \n\) QM = ^ / d 3 pu(p, \')a 03l5 u(P, A)|<Kp)| 2 • (70) 

Here, |nA) is a neutron in the S z — A state, the <p(p) are used in the wave 
packet description of the momentum eigenstates as discussed in Appendix 
A, and the u(p, A) are the neutron spinors. Choosing A = A' = 1/2 we have 



(nl/2\ f d 3 x 3 J eff \nl/2) QM = ^ . (71) 



dr, 

e 

We evaluate the two matrix elements needed and find 



t „ v .8 0.8623 

{n\uj 5 udj d\n) = i- ^ R2 (72) 

/ i ou 7 n x 4 0.8623 

(n\ua°"ud l5 ^d\n) = --^ . (73) 

Using these matrix elements, we obtain the following expression for the neu- 
tron EDM: 

A TVPC \A 2 TVPC RIJ 32™^ ° g \Ml) {l > 

x^[g d v(g u A + g d A) + ^g d Agv] (75) 

Since g\ + g d A ~ in the Standard ModelQ, this contribution to the neutron 
electric dipole moment is approximately proportional to g^gy- 

The second way TVPV operators contribute to d n is by mixing states of 
opposite parity into the neutron ground state (Fig. 4a). The lightest state 
which may contribute is the Nn S-wave. In the chiral limit, its contribution is 

6 The sum is exactly zero at tree-level. 
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dominated by the loop diagrams of Fig. 8(b). The TVPV NNtt couplings are 
just the g(p' ', generated from the TVPV quark operators as in Fig. 8(a). The 
g^' also contribute to the atomic EDM via the meson-exchange interaction 
of Fig. 1(c). In what follows, we relate the g^' to the purely hadronic parts 
of the operators in Eq. fl45|) . 

Following Ref. [[Uj, we carry out the calculation in the factorization 
approximation while using the on-shell form of the TVPV operators given 
in Eq. (f46"l) . First, we consider the TVPV NNtt coupling. The only Dirac 
structures which give rise to pion- vacuum matrix elements are 75, and 757^. 
Therefore the structure in the third line does not contribute. The last two 
structures in Eq. (f46"D do not contribute in the factorization approximation. 
This conclusion follows from symmetry arguments, p^p^a^ = 0, and the 
equations of motion p^N^^N = 0. For the remaining structures the PCAC 
relation 

-(0|m7 m 7 5 m - ^^dlTr ) = if^ exp(-ip • x). (76) 

Taking the divergence of Eq. ([F|) and using the consequences of isospin 
symmetry 

(0|n 75 u + d l5 d\ tt°) = (77) 
yields the following two matrix elements needed in the calculation: 

(0|m7 5 -u|7t ) = -i exp(-ip ■ x) (78) 

m u + m d 

(0| d lb d\ix°) = % m2?rfn exp(-ip ■ x) . (79) 
m u + m d 

In addition, we require the nucleon matrix elements of the light quark 
scalar densities. From the quark model and ttN sigma term, one obtains 

(n\uu\n) = (p\dd\p) ~ 5 (80) 

(n\dd\n) = (p\uu\p) ~ 6 (81) 

Using these results, we obtain for the TVPV neutral pion nucleon couplings 
An7c°n and Bpir p, the strength of the couplings are: 



A = lr a7 ""' oo,r 2 .2 tog ( T72 ) :— fV 36 9v9a + ^9 u v (g u A + 9 d A ) 



3 Al VPC 32ns 2 v c 2 v \Ml m u + m d 



10 



36^ 

(82) 



26 



U m l a , ( I** \ m d 



2 



3 A^ yPC , 32vr,s^c% \M 2 Z ) m u + m d L 5 

(83) 

Here again, since + gr^ = at tree level in the Standard Model, the terms 
proportional to gyg\ make the largest contribution. 

In order to determine the couplings of the charged pion with the nucleon, 
we Fierz transform Eq. (|46D . The result is listed in Appendix B. As discussed 
in Appndix B, we estimate the strength of the Cpir~n coupling using the first 
two terms in the Fierz transformed expression. The remaining terms do not 
contribute in the factorization approximation. In this case, we require the 
following matrix elements: 

f 2 

(0\u l5 d\n-) = (n-\d l5 u\0) = 1V2 U ^ ; (84) 

m u + m d 

(n\du\p) = (p\ud\n) rs I. (85) 
The coefficient can then be written as, 

3 \A 3 TVPC ) 32nsWw g Wz) 

md ( ' -u d 1 u u 1 3 d u 



m ,, + ,n d y^-2^ + 4^) < 86) 

The last piece 3gyg^/4 comes from applying equations of motion to the 
derivative terms, as discussed in Appendix B. 

Expressed in terms of A, B and C, the PVTV pion- nucleon couplings 
from Eq. (|3]) are 

gW = ( B -A)/Q + V2C/3 (87) 
g^ = {A + B)/2 (88) 

gW =(B- A)/6 - ^= (89) 

In the chiral expansion of the EDM, the dominant contribution to d n 
from the TVPV ttNN interaction arises from the loops in Fig. 8(b), where 
the intermediate state contains a p7i~. In this case, only the constant C 
contributes. The loop calculation, first performed in Ref. JfCj, yields 



d n _ C g-n-NN itin (90) 
e \/2~4:7r 2 mN rn n 
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Written explicitly in terms of the TVPC scale, Atvpc, the neutron electric 
dipole moment takes the form, 




a 



3 32ns 2 w c 2 w 4vr 2 



log- 



m N 



±9v9a - \dv9 U A + -/y9a 



(91) 

where we have kept only the leading, non-analytic loop contribution pro- 
portional to log(m n / mjsr) and where g n NN is the strong 7riViV coupling. We 
note that an evaluation of the analytic contributions has been performed us- 
ing sidewise dispersion relations in Ref . [ 31 1 . In addition, loop contributions 
involving the nn intermediate state have been considered in Refs. |32], |33 



D. Tree-level contributions 

There exists one way in which the d = 7 TVPC operators may contribute 
to d n without the consideration of loop effects. As shown in Fig. 9, the 
operator 0^2 generates a tree-level contribution when the Z° is exchanged 
and couples to the second quark's axial vector current. Naively, one might 
expect this contribution to compete with the one-loop quark EDM generated 
by 0^2 m Fig- 2. The calculation of the process in Fig. 9 is tedious but 
straightforward. Using the quark model, we obtain 



d n 
e 



1 



A 7 



TVPC 
2.04 

2 



2s w 



My, 



An 



0.34 



(92) 



X 



E 



a; 



lo(x) 2 fo( x )ji( x ) - jo(x)ji(x) 2 j[(x) - j (x)ji(xf/x 
(nl/2|u + (m', a')u(m', a)cft(m, f3')d(m, p)\nl/2) , 



m,m' ,a, 



where the sum runs over all spin (rn,m') and color (a,/3, . . .) indices. Eval- 
uating the integrals and quark model contractions yields 



dn 
e 



A 



TVPC 



W 9 A 
2s w 



M, 



A 7 



0.162 



MlM w Ri 



(93) 



Comparing the result in Eq. ( PB| ) with one-loop amplitude in Eq. 
we observe that the tree-level contribution contains the suppression factor 

0.162 



8.64 x 10 



-12 



(94) 
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The difference results from the contributions of high- momentum (p ~ M z ) 
intermediate states in the loops of Fig. 2. The tree-level process, in con- 
trast, is dominated by states with momenta p^A QCD ~ 0.002M Z . Since 
the TVPC operator is has dimension seven, and since the momentum trans- 
ferred through Z°-boson propagator in Fig. 9 is ^A QCD , one would expect 
a suppression factor of (A QCD /M Z ) 4 ~ 1CT 11 for the process of Fig. 9. 



6 Experimental Constraints 

The foregoing analysis allows us to make a quantitative connection between 
various observables and effective coupling constants on the one side and the 
EFT of TVPC interactions of quarks and gauge bosons on the other. In this 
section, we use those relationships to derive bounds on the TVPC interac- 
tions. It is also instructive to interpret these bounds in terms of the mass 



scale, A TVPC . To that end, we adopt the parameterization of Eqs. ([34] , |35| ) 
and present constraints in terms of A TVPC and k. 

In what follows, we use experimental d e and d n limits directly as well as 
the limits on g p and the g^f 1 derived from several sources. The corresponding 
limits on these quantities are given in Tables 1 and 2. It is also useful to 
convert our expressions for g p and the g^ into numerical form. In the case 
of g p we have 

g p « -Ct {M z /A TVPC f x (2 x 1(T 10 ) . (95) 

Similarly, the combinations of the g^ relevant to d n and c? A ( 199 Hg) are, 
respectively, 

# - # « -C 7a {M z /A TVPC f x (6 x 10- 12 ) (96) 

for d n and 

# + g<» + 2# « -Ct {M z /A TVPC f x (3 x 10~ n ) (97) 
for d A ( 199 Hg). 

Scenario A. Under this scenario, the first two terms in Eq. (|T]) vanish, so 
that the leading terms in the expansion of the EDM are the 0(1/A^, VPC ) con- 
tributions from the TVPC operators. As in the case of the effective hadronic 
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couplings, it is useful to express our relationships between the EDM's and 
the TVPC interactions in numerical form. 

For the elementary fermion EDM's arising from the one-loop graphs con- 
taining O^c (Fig. 2), we require a choice for /i. Since the typical momentum 
of a quark in the neutron is ~ A QCD , we take \x = A QCD . The appropriate 
choice for d e is probably smaller, on the order of the typical momentum of 
an electron bound in a heavy atom. Since the /i-dependence of df is only 
logarithmic, however, the precise choice of scale is not decisive. To be con- 
servative, we also use \i = A QCD for d e . The corresponding results are 



d e 



and 



{M z /A TVPC f x (4 x 10" 17 cm) (98) 

C{{ {^a/s w ){M z /A TVPC f x (4 x 1(T 17 cm) 
C//' {M z /A TVPC f x (8 x 10~ 18 cm) 

(99) 



- - {\A ~ k) x (0-88) x *± (100) 
e 3 3 e 

i z ~i C\t I A \% ^ i a v in -i7 



Cf c 7 (M Z /A TVPC ) J x (6 x 10- 17 cm) 
C//' (M Z /A TVPC ) 3 x (1 x 10 



-17 N 

cm) , 

where we have used the naturalness relation of Eq. ( p5|) to express Cj^ in 
terms of C(l and where the 0.88 factor in Eq. ( |100| ) is the value of the bag 
model integral in Eqs. (^) and ([5(]). Since the two- loop effects involving 
Ojl are suppressed relative to the one-loop 0^2 contributions, we do not 
give explicit numerical formulae for the former. 

For the contribution to d n from the TVPV 7-four quark interaction of 
Fig. 4b, generated by the loops of Fig. 5b, we find 

^ ~ Cf a (M z /A TVPC f x (5 x 10- 25 cm) , (101) 

while the process of Fig. 9 yields 

d n 



rsj 

e 



C£f {M z /A TVPC f x (4 x lO" 28 cm) , (102) 
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where we have again used Eq. (|35|). 

We now apply these expressions to the experimental limits given in Tables 
1 and 2. The results are listed in Table 3. Clearly, the one-loop elementary 
fermion EDM's generated by radiative corrections to OjJ yield the most 
stringent bounds on A TVPC . Whether one uses d n or d e , the lower bounds 
generated by this mechanism are roughly three orders of magnitude larger 
than those obtained from the many-quark effects. The reason for this dif- 
ference is clear. The one-loop graphs are dominated by high-momentum 
(p ~ M z ) intermediate states, whereas the many-quark matrix elements are 
governed by physics at scales p ~ A QCD . Since the dimension of the TVPC 
operators is greater than that of the EDM by two, the EDM must contain at 
least two powers of the relevant mass scale. Hence, the one-loop (9^7 effects 
should be at least [M z j A QCD ) 2 ~ 2 x 10 5 larger than the many-quark effects 
in d n . 

We emphasize the comparison between the one-loop elementary fermion 
EDM limits and those obtained from d A ( 199 Hg). In the latter case, the 
strongest bounds are derived from the combination of the appearing 
in Eq. (|97|) . Since these couplings scale as 1/A 3 VPC , one would require 
an improvement of nine orders of magnitude in the experimental limits on 
<i A ( 199 Hg) in order for the atomic EDM to compete with d e and d n f\ We 
also note that were it not for the one-loop elementary fermion limits, the 



diagrams of Fig. 5b - required by gauge invariance but omitted in Ref. [[16 
- would yield the most stringent constraints on new TVPC interactions un- 
der scenario (A). By far the weakest bounds follow from the extraction of g p 
from <i A ( 199 Hg). 

The limits in Table 3 can be translated into an expectation for ax under 



Scenario (A). Using the general scaling arguments of Ref. (T8| which imply 

, 3 „„/ / M z \ 3 / p \ 3 

^■TVPC f V M z , 

and using the naturalness assumption of Eq. (|35|) we have 

3 

-16n 



a T ~ (7 x 10- 18 ) x (j^J • (104) 



7 An improvement of this magnitude, however, would likely render the atomic EDM as 
the most precise probe of d n . 
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For a low energy TVPC process with p ~ 1 GeV, the Scenario (A) EDM 
limits of Table 3 imply that the size of the effect should be roughly 10~ 15 . 
Current direct search limits are at the 10 -3 level for a-r- 

Scenario B. As argued above, experimental EDM limits do not apply to 
TVPC interactions when PV persists at short distances. In this case, one 
must rely on direct TVPC tests, from which limits on g p may be extracted. 
The most precise determinations of g p from TVPC observables are obtained 
from charge symmetry breaking (CSB) terms in np scattering cross sections 
[ PI and fivefold correlations (FC) in the transmission of polarized neutrons 



through spin aligned Holmium |T0|. The limits on g p obtained from these 
observables are given in Table 2. In addition, a proton-deuteron transmission 
experiment proposed for COSY would be sensitive to g p at the 10~ 3 level 
Further improvements in the np CSB limits may also be possible |14 . 



The bounds on A TVPC from direct TVPC searches are given in Table 4. 
At present, the data imply new TVPC physics could arise at scales as light 
as a few GeV. Ideally, future experiments would push these bounds closer to 
the weak scale. Such a sensitivity would be comparable to the present and 
anticipated new physics sensitivities of atomic PV and PV electron scattering 
experiments Q. Achieving such sensitivity, however, would be a daunting 
task, requiring at least six orders of magnitude improvement in precision 
beyond the present state of the art. 



7 Conclusions 

Measurements of EDM's provide one of the most powerful probes of possible 
new physics. Thus far, null results for d e , d n , and d A lead to tight constraints 
on a variety of new physics scenarios. In this paper, we have studied the 
implications of EDM's for TVPC new physics. We have also developed re- 
lationships between effective TVPC and TVPV hadronic couplings, on the 
one hand, and TVPC interactions involving quarks and gauge bosons, on the 
other. This relationship has not been systematically delineated in the past. 
Consequently, some confusion about the respective implications of EDM's 
and the values of hadronic couplings extracted from direct TVPC searches 
has ensued. We believe that our analysis has helped clarify these implica- 
tions. 
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While the possible origin of TVPC interactions in the context of a renor- 
malizable gauge theory has yet to be delineated, one may, nevertheless, ad- 
dress the issue using the framework of EFT. This framework affords a sys- 
tematic method for treating nonrenormalizable interactions and for carrying 
out the associated phenomenology. Since the operators we have treated here 
have d > 4, EFT is the natural framework for performing our analysis. As 
emphasized above, a key ingredient in the application of EFT to this prob- 
lem is the maintenance of a scale separation. The short- distance (p > A TVPC ) 
TVPC physics about which we are ignorant is parameterized by the a priori 
unknown coefficients of the nonrenormalizable effective interactions. Long 
distance (p « A TVPC ) physics may be treated explicitly in the guise of 
loops and many-body matrix elem ents involving degrees of freedom having 
masses and momenta below A TVPC . As we illustrated for the example of 
the cut-off regulator used in the loop calculations of Refs. WQ, [T7J], a failure 
to adhere to this scale separation can lead to disastrous results. Indeed, in 
that example, the blurring of scale separation implies (erroneously) that the 
EDM is proportional to an infinite series of unknown coefficients of nonrenor- 
malizable operators of arbitrarily high dimension. Without the truncation 
scheme implied by the EFT scale separation, one is unable to learn anything 
from experiment about a given TVPC interaction (or even finite set of such 
interactions). 

As argued at the outset of our paper, the EFT scale separation implies 
one must consider the low-energy consequences of TVPC new physics un- 
der two scenarios: parity-restoration below A TVPC (Scenario A) and parity- 
restoration above A TVPC (Scenario B). In the case of Scenario A, the EDM is 
dominated by PV radiative corrections to d — 7 TVPC operators. The most 
significant impact arises from the EDM of elementary fermions, generated 
by a one-loop graph involving 0^2 ■ Interpreted in terms of a mass scale, 
this mechanism, taken with the experimental limits on d e and d n , imply that 
A TVPC > 100-250 TeV when the TVPC new physics is "strong" (k ~ 1). 
The impact of many-quark effects - which contribute both directly to d n as 
well as to both d n and <i A ( 199 Hg) via the couplings g p and g^>' - is consider- 
ably smaller than that of the one-loop elementary fermion EDM. Indeed, for 
the many-quark mechanism to compete with that of the elementary fermion 
EDM, one would require an improvement in the limits on <i A ( 199 Hg) of nine 
orders of magnitude. The impact of direct search limits is even weaker under 
this scenario. We conclude, then, that d e and d n provide the most power- 
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ful probes of new TVPC physics when parity symmetry is restored at short 
distances, while the atomic EDM measurements are better suited to con- 
straining other new physics scenarios. 

Under Scenario B, the implications of EDM's are ambiguous at best. The 
existence of parity- violation at short distances implies that d < 7 TVPV op- 
erators contribute to the EDM along with the radiatively corrected d > 7 
TVPC operators. Consequently, no single TVPV observable constrains the 
latter unless one invokes additional assumptions. In contrast, the relation- 
ship between TVPC observables - such as CSB np scattering cross sections 
or FC in neutron transmission - is not clouded by the additional unknown 
constants which enter the EDM. In this case, however, the limits obtained 
from direct TVPC searches are rather weak: A TVPC > 1 GeV (for k ~ 1). 
As a long term benchmark, one would ideally search for new TVPC physics 
at least up to the weak scale. Achieving this goal would require an improve- 
ment in the precision of direct TVPC searches by six orders of magnitude. 
While achieving such improvement would seem formidable, any intermediate 
progress would constitute a welcome step. 
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A Quark model matrix elements 



Here, we evaluate the matrix elements appearing in Eq. (|57|) . We define the 
following quantities: P = pn + Pp, q = P P = Pn — Pp, e is the p polarisation. 
Since the TVPC pNN Lagrangian is linear in the p momentum operator, 
we retain only those terms in the matrix element products linear in q. The 
fourth matrix element product in Eq. ( |5?j ) is quadratic in q since the nucleon 
matrix element vanishes to first order in q. Thus, the last producet does not 
contribute to g p . 

For the remaining matrix element products in Eq. ([57]), we first require 
the p to vacuum matrix element, parameterized in the standard way as 



2 

-m; 



(0|«7 M d|p-> = v^yV^ (105) 
Jp 

where /p/47r w 2.5. The rho matrix element with the derivative may be 
evaluated using wave packets in the quark model, using the approach of 
Donoghue and Johnson p5] . We begin by assuming the following structure 
for the matrix element, 

(0|u^" M d|p-) = (ae„ + bp^)e ipx (106) 

where p^ is the momentum of the rho-meson and is the polarization vector. 
We solve for the coefficients, a and b. The quark model matrix element may 
be written using wave packets as 



d 3 p' 

Multiplying each side of the equation by / <i 3 xexp(— ip ■ x)/(2ti) 3 , results 

in 

d 3 x ._, -^ /n |_^ _<j>{p), 

(2?r) ..... ^ - _ p 

The expression for <p(p) may be obtained from the vector current matrix 
element Eq. ( |105| ). 



(0\Ud,d\p-) QM = J ^{p'^ud ,d\p-) (107) 



f d 3 x < (t)(v) 

J (2ttP %)(Q\ud ^d\p~) QM = ^—(ae^ + bp^). (108) 



— — -exp(-ip- x){0\unf»d\p-) QM = (f)(p) / P r- (109) 
( 27r ) fpOOpVZ 
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We use the /i = 3 component of the above equation to find an expression for 
<f)(p). This, taken together with the vector current matrix element produces, 



(110) 



Sxe-^ s ^\u lz d\p-) QM 



Upon evaluating both quark model matrix elements, we find 6 = 0, while, 



m 2 2 
X3 



du dl 
—I — 2ul r — u— 
dr dr 



r 2 dr/ / (u 2 - l 2 /3)r 2 dr. 



(Ill) 



Here, u and I are the upper and lower components of the Quark Model wave 
functions, respectively. The final expression is 



\ud p d\p ) 



1.05 r-m 2 p . _ 



■ip-x 



(112) 



In this expression, R p = 3.3 — 3.5 GeV -1 is the radius of the rho meson in the 
bag model. In the following we shall also use R n = 1 fm for the bag model 
radius of the nucleon. The nucleon matrix elements also may be evaluated 
with the quark model. 



(n\d-i^u\p) = u n (^ + 0.l68iR n a^q u )u p . 



(113) 
(114) 



{n\d d u u\p) = O.Q3Aiuju n / y I/ Up 

where uj = 2.04/ R n is quark energy inside the bag. 

For the third matrix element product, we obtain using similar methods 
and Eqs. ((48], (49]) 

^(0\u(d^u - *d v 7n)d\p~) = \ 5 + +3g A r^/2f ~ e A*5") 
(n\da^u\p) = ^(1 + 9 -^)u N a^u P + ■■■ , 



where the + • ■ ■ denote terms higher order in q. 
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B Fierz transformed four-quark operators 



In calculating the TVPV 7r ± -nucleon couplings, we require the Fierz trans- 
formed form of the amplitude in Eq. ([461) . The transformed amplitude is 



J^FIERZ _ -in 



a • - 2 



^■TVPC ^ 7rs W c W ^ 1V1 Z' 

im r (Ag{g f v - -g( r 9 f A )(U f jBUfU r U f + U f U r U fl5 U f ) 

+im f ,(2 g y v + l -g f v g f A ) (U fll ,U f U f ^ l5 U f - U fl ^ 5 U r U f ^U 

-\ m f'9v(.9A + gAWfcr^UfUf^a^Uf + ^m f/ g{ / g f A Ufy 5 a^U f U f a lliy Uf 

+i(p' f +p f r(2g A 'g{; - g£g A + \g f A gC)\{UfU f U f H^U f - U flfl y 5 U f ,U f U f ,) 

7 1 

+i(p' f + VsY^9a9v ~ 9v9a + ^9A9vh{U flh U fl U n ^Jf + U n „U f U n5 U f ) 



2 A v M 

(p' f + p f T{2 g y; - g£g f A + L{ g £)hu f a^ l5 U f ,U f ,YU f - U n »U r U r a^U f ) 



-{p' f + p f y{2 9 y; - g£g f A + lg A g f v)~ A {Un V l*U fl U r a, u U f + U f a^UfU f ,Y l5 U f ) 
3 



2 4 

+^9 f A9v i (p'f +Pf') f *(Un»U f U f rj 5 U f + tJ f j 5 U f U f ^U f ) 



8^ 
3 ' 

--^gigvKp'f + pMUfy^UfiUfiUf - u f u fl u f ,^ 5 u f ) 

-ggigvip'f'+Pf^iUfcr^UfUf.j^Uf + U f j v j 5 U f >U f a^U f ) 

+^9A9v(p'f' + Pf')^Ufj v Uf>U r j 6 <T vlt Uf - UfYsa^UyUf^Uf)} (115) 

Only the first two terms on the RHS of Eq. ( |115| ) contribute in the factor- 
ization approximation. The terms on the second line, such as Ufj^UfiUfi^^U f 
will not contribute in this approximation. The reason is that the pion-to- 
vacuum matrix element will be proportional to the pion momentum, q, which 
equals the change of the momentum of the nucleon. When contracted with 
the nucleon matrix element of the vector current, the result vanishes by cur- 
rent conservation (neglecting the small n-p mass difference). The terms on 
the third line will not contribute for symmetry reasons. 
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The remaining terms involve derivatives on the quark fields. In each case 
the derivative pair is split, so one derivative acts on a quark field in one quark 
bilinear and the second derivative acts on the second. Some of these terms 
may be rewritten using equations of motion for the on shell quarks. Doing 
so yields new contributions to the first four lines of the RHS of Eq. ( |115| ), 
requiring an adjustment in the coefficients. The remaining terms involve 
derivative operators, such as 

Ufp'^UfUfj^Uf and U f p' f ^ 5 U f U f ^U f . (116) 

Factorization matrix elements of these operators vanish by either symmetry 
or current conservation. 

Finally, tensor structures such as U^p'f^UfUfa^j^Uf cannot produce 
a pion- vacuum matrix element, while those such as Lffp'^j^Uf/Ufia^Uf 
will vanish from symmetry considerations. 
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Table 1: Present experimental limits on electric dipole moments. 
Observable Experimental limit (e-cm) Ref. 



d e 


< 4 x io~ 27 


U 






< 8 x 1(T 26 


29 




^( 199 Hg) 


< 1.3 x 10" 27 







Table 2: Experimental limits on TVPC and TVPV hadronic couplings. The 



value of h°§% is the "best value give in Ref. 
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Coupling 



Experimental limit 



Observable 



9 P 
9 P 



#1 



S( 2 )| 



< 9.3 x 10~ 3 

< 0.53 x lQ-\h^ H N /K NN ) 

< 5.8 x 1(T 2 

< 6.7 x 10~ 3 

< 5.7 x 10~ 12 

< 1.8 x KT 11 



ci A ( 199 Hg) & atomic PV 
d n 

FC in neutron transmission 
CSB in np scattering 
d n 

^( 199 Hg) 



Table 3: Limits on new TVPC interactions derived from EDM's under Sce- 
nario (A). 



Observable Lower bound on A TVPC / k 2 ^ 3 (TeV) 



Mechanism 



single electron 0^2 1°°P 
single quark loop 
Figs. 4b, 5b 
7r-loop, Fig. 8b 

Fig. 9 
# (Fig. lc) 
g p , atomic PV (Fig. lb) 





d e 


260 




dn 


110 




d n 


0.39 




d n 


0.21 




dn 


0.036 


d A 


: i99 H g ) 


0.25 


d A 


: i99 H g ) 


0.0006 
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Table 4: Limits on new TVPC interactions derived from direct searches under 
Scenario (B). 

Coupling Lower bound on A TVPC / k 2 ^ (TeV) Observable 

g p 3 x 1CT 4 CSB in np scattering 

g p 7 x 1CT 4 FC in neutron transmission 
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Figure 1: Representative contributions to atomic EDM: (a) TVPV nuclear 
effect involving TVPC and PV pNN interactions; (b) TVPC nuclear effect 
plus atomic PV; (c) long-range TVPV nuclear effect involving TVPV ttNN 
coupling. Open circle denotes strong meson-nucleon coupling; crossed circle 
gives TVPC coupling; open square is TCPV coupling; and crossed square is 
TVPV coupling. 

Figure 2: One loop contribution to elementary fermion EDM from the 
TVPC operator 0? z . Coupling symbols are as in Fig. 1. 

Figure 3: Two loop contributions to the elementary fermion EDM involving 
the TVPC operator O f / . Symbols are as in Fig. 1. 

Figure 4: Contributions from four-quark TVPV operators to d n : (a) second 
order contribution involving a mixture of opposite parity states \n) into neu- 
tron; (b) first-order contribution arising from 7-four quark TVPV operators. 
Symbols are as in Fig. 1. 

Figure 5: PV weak radiative corrections to 7a , generating d = 7 TVPV 
operators. Symbols are as in Fig. 1. 

Figure 6: Contribution to d n from individual quark EDM's. Symbols are as 
in Fig. 1. 

Figure 7: (a) Contributions to g p from Oj[ ; (b) contribution to d n arising 
from TVPC pNN and PV ttNN interactions. Symbols are as in Fig. I. 

Figure 8: (a) Contributions to gW' arising from TVPV many-quark opera- 
tors; (b) leading-oder contribution (in m n ) to d n arising from TVPV 7iNN 
interaction. Symbols are as in Fig. 1. 

Figure 9: Tree-level, many-quark contribution to d n generated by ■ 
Symbols are as in Fig. I. 



44 



+ 




+ 




+ 





Q. 



